The QCD vacuum as a disordered medium: 
A simplified model for the QCD Dirac operator 
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We model the QCD Dirac operator as a power-law random banded matrix (RBM) with the ap- 
propriate chiral symmetry. Our motivation is the form of the Dirac operator in a basis of instantonic 
zero modes with a corresponding gauge background of instantons. We compare the spectral corre- 
lations of this model to those of an instanton liquid model (ILM) and find agreement well beyond 
the Thouless energy. In the bulk of the spectrum the (dimensionless) Thouless energy of the RBM 
scales with the square root of system size in agreement with the ILM and chiral perturbation theory. 
Near the origin the scaling of the (dimensionless) Thouless energy in the RBM remains the same as 
in the bulk which agrees with chiral perturbation theory but not with the ILM. Finally we discuss 
how this RBM should be modified in order to describe the spectral correlations of the QCD Dirac 
operator at the finite temperature chiral restoration transition. 
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In recent years the relation between spontaneous chiral 
symmetry breaking (S^SB) in QCD and the phenomenon 
of conductivity in a disordered medium has been investi- 
gated in the literature UJ- L3- Lj| • In the latter, conductivity 
is produced by electrons that although initially bound to 
impurities may get delocalized by orbital overlapping. In 
QCD the quark zero modes initially bound to an instan- 
ton get delocalized due to the strong overlap with other 
would-be zero modes and consequently chiral symmetry 
is broken. In the case of atoms the overlap is effective 
only among nearest neighbors due to the exponential de- 
cay of the electron wavefunction. However, in the QCD 
vacuum, the decay is power-law and long range hopping 
is possible. Thus even if one assumes a perfect trade be- 
tween impurities and instantons and between electrons 
and quarks, the associated Anderson model should posses 
long range hopping. Such models have already been dis- 
cussed in the literature 0, 0. The main conclusion of 
these works was that power-law hopping may induce a 
metal insulator transition even in one dimensional sys- 
tems if the exponent of the hopping decay matched the 
dimension of the space. Similar findings were also re- 
ported in a related model, a random banded matrix with 
a power-law decay 0, Q] . 

Random matrix techniques have already been used in 
the context of QCD. In the infrared limit the eigenvalue 
correlations of the QCD Dirac operator do not depend on 
the dynamical details of the QCD Lagrangian but only on 
the global symmetries of the QCD partition function || . 
Thus random matrices with the correct chiral symmetry 
of QCD (termed chiral random matrices) accurately 
describe the spectral properties of the QCD Dirac oper- 
ator up to some scale known as the Thouless energy. For 
larger energy differences dynamic features become im- 



portant and the standard random matrix model ceases 
to be applicable. 

In this letter we study an improved matrix model for 
QCD which reproduces features of the eigenvalue spec- 
trum beyond the Thouless energy. We want to incorpo- 
rate the phenomenological fact that the matrix elements 
of the QCD Dirac operator in a basis of zero modes decay 
as a power of the instanton-anti-instanton distance into 
the random matrix approximation. Wc thus propose a 
chiral RBM as an effective model of the QCD Dirac op- 
erator. We first give a brief introduction on the role of 
instantons in the QCD vacuum. Next we introduce the 
chiral RBM and describe its main properties at both the 
origin and bulk of the spectrum. In the latter case an 
analytical treatment is possible. Then we compare the 
spectral correlations of this model with the results of an 
instanton liquid simulation. Finally we discuss how the 
RBM can be modified to describe the spectral correla- 
tions of the QCD Dirac operator at the finite tempera- 
ture chiral restoration transition. 



SxSB FROM INSTANTONS 

The discovery of instantons has had a large impact on 
our understanding of non-perturbative aspects of QCD 
[lfl | (for a modern review see 0). Instantons are clas- 
sical solutions of the Yang-Mills equations of motion in 
Euclidean space. An important property is that the Eu- 
clidean QCD Dirac operator has an exact zero eigenvalue 
in the field of an instanton. The spectral properties of the 
low lying modes of the Dirac operator are thus controlled 
by these non-perturbative configurations. 

Unfortunately the construction of a consistent QCD 



2 



vacuum based on instantons faces serious technical diffi- 
culties. Exact analytical multi-instanton configurations 
are hard to obtain since the Yang-Mills equations of mo- 
tion for QCD are nonlinear and therefore a superposition 
of single instanton contributions is not itself a solution. 
Additionally quantum corrections may spoil the semi- 
classical picture implicitly assumed of a QCD vacuum 
composed of instantons well separated and weakly inter- 
acting. These problems have been overcome either by 
invoking variational principles \l2\ or by phenomenologi- 
cally fixing certain parameters of the instanton ensemble. 
The latter case, usually referred to as the instanton liquid 
model (ILM) [13J, yields accurate estimates of vacuum 
condensates and hadronic correlation functions just 
by setting the density of instantons to be N/V « 1 fm -4 
and the mean size p w 1/3 fm. Lattice simulations have 
also supported the picture of a QCD vacuum dominated 
by instantons 15] . 

An advantage of these models is that they satisfactorily 
explain how chiral symmetry is spontaneously broken in 
the QCD vacuum |l2j. The chiral condensate (the order 
parameter signaling SxSB) (ipijj} is related to the spec- 
tral density p(e) of the QCD Dirac operator around zero 
through the Banks-Casher relation 
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where V is the space-time volume. The ILM provides a 
phenomenological model for these low energy modes of 
the QCD Dirac operator. In a basis of N/2 instantonic 
zero modes ipi(zi) and N/2 anti-instantonic zero modes 
ipA(zj) the matrix elements of the Dirac operator take 
the form (for zero quark mass) 
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with the N/2 x N/2 overlap matrix given by 

Tia = (lpl(Zi)\ft\lpA(Zj)) ~ piPj/\Rij\ c 
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Zj between the center 



for large separations Rij = z 
of instanton i and anti-instanton j with pi and pj their 
sizes. Due to the chirality of the zero modes only matrix 
elements connecting an instanton with an anti-instanton 
do not vanish. Physically the amplitude of the matrix el- 
ements Tja represents the probability for a quark to hop 
between an instanton and an anti-instanton. An isolated 
(anti-)instanton would cause an exact zero mode but this 
degeneracy is lifted through overlap with neighboring in- 
stantons. Thus the would-be zero modes are effectively 
split around zero. As more (anti-)instantons are added a 
continuous band spectrum is formed with a spectral den- 
sity finite at zero. As mentioned previously, it is precisely 
this dynamically generated non-zero spectral density that 
causes, through the Bank-Casher relation, a finite value 
of the chiral condensate. 



THE CHIRAL RANDOM BANDED MODEL 

Here we study the spectral properties of an ensemble 
of chiral random Hermitian N x N matrices given by 
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where C is a N/2 x N/2 complex matrix with indepen- 
dently distributed Gaussian variables with zero mean. 
The variance of the matrix elements CL- are chosen to de- 
cay as a power of r — which measures the distance 
from the diagonal. Since the ILM uses periodic boundary 
conditions we use a periodic form of the power-law decay 
[T?! given by 
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where a and b are real parameters. The choice of com- 
plex matrix elements corresponds to a matrix model with 
a unitary symmetry which is appropriate for QCD with 
the phenomenologically relevant ST/ (3) color group. Due 
to the chiral symmetry, the eigenvalues of J3J come in 
pairs of ±€i. This feature induces an additional level re- 
pulsion around zero which results in different spectral 
correlations for eigenvalues near zero (the origin) and 
away from zero (the bulk). 

In the bulk the spectral correlations should not be af- 
fected by the block structure and should coincide with 
the non-chiral version of © which has been intensively 
studied in recent years [3, Il7| . The use of the supersym- 
metry method |18| permits an analytical evaluation of 
both spectral properties and eigenfunction statistics Q 
in a certain region of parameters. In the thermodynamic 
limit the eigenfunctions are multifractal for a = 1 and 
localized (delocalized) for a > 1 (a < 1) respectively 0- 
The spectral correlations in the g = E c /A >> 1 (E c is 
the Thouless energy and A is the mean level spacing) 
limit can be expressed through the spectral determinant 
of a classical diffusion operator 0] . The two point corre- 
lation function is defined as i£a(s) = A 2 (p(e)p(e+sA)) — 1 
where p(e) is the density of states at energy e and the 
average is over an ensemble of RBM. For the unitary en- 
semble (our case) 

„ , s 1 d 2 , D(s,q) cos(2ns) . 

R ^) = -^g^^^ M1 + ^Z^ i D(s,g). (6) 
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Due to the power-law decay, the spectral determinant 

D(s,g) = IW 1 + ^hl)- 1 ( = sM^lcorre- 
sponds with a process of anomalous diffusion y}. For 
1/2 < a < 1 the dimcnsionless conductance increases 
with the system size as g = C a (b)N 2 ~ 2a with C a (b) a 
known constant. The scaling of g thus resembles that of 
a weakly disordered conductor in d = 2/ (2a — 1) dimen- 
sions [7|. 
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RBM analytic g = 21 

RBM N = 640 
RBM analytic g = 30 
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FIG. 1: Spectral rigidity for the RBM in the bulk. Points 
correspond to the numerical simulations of the chiral RBM 
for 6 = 1. Lines correspond to the analytical result of the 
non-chiral RBM © for the given values of g. 




FIG. 2: Spectral rigidity for the ILM and RBM in the bulk. 
Points correspond to numerical simulation of the ILM. Lines 
correspond to the analytical result of the non- chiral RBM © 
for the given values of g. 



The above results were derived for the non-chiral RBM. 
As mentioned above they are also expected to describe 
the chiral RBM J3J in the bulk region. We have chosen 
the scale of b in (J3J to agree with results for the non- 
chiral model. Strictly speaking, the factor C a (b) was 
derived for b ^> 1 and it is not clear how the result 
should be modified for smaller b. Below we will show 
that in fact the analytical results agree very well even 
at b = 1. Unlike in the bulk, analytical results for the 
spectral correlations close to the origin are not known. 
They can be obtained by modifying the supersymmetry 
method to account for the chiral structure. This issue 
will be postponed to a future publication |20| and in this 
letter we will rely on numerical calculations for the study 
of the spectral correlations close to the origin. Lastly we 
remark that though the power-law decay in the instanton 
liquid is not in principle related to the matrix index, as in 
our matrix model, numerical simulations show that both 
models yield similar results We stick to Q due to 
the availability of analytical results. 

Given the known results for the non-chiral RBM we 
can now choose the parameter a. Recall that our moti- 
vation to propose a random banded model is the decay of 
the ILM overlap matrix elements J3J . It was shown in [j| 
that the spectral properties of systems with power-law 
hopping are similar in different dimensions provided that 
the decay exponent equaled the dimension. Since in the 
ILM the decay exponent (three) is less than the dimen- 
sion (four) we expect this to map onto a ID RBM model 
with a < 1. We choose a = 3/4 because the volume de- 
pendence of the dimensionless conductance (also called 
dimensionless Thouless energy) in this case g s» 1.17 'VbN 
coincides with what is expected for QCD according to 
chiral perturbation theory |2]J. As mentioned above, the 
spectral properties of the RBM at a = 3/4 are similar to 
those of a disordered conductor in four dimensions. 



RESULTS 

We generated sets of matrices for the RBM and the 
ILM and calculated their eigenvalues in order to compare 
the spectral correlations. The results of each simulation 
are averaged over 10 4 configurations. For both models 
the simulations were done in the quenched approximation 
which allowed us to use values of N ranging from 320 up 
to 1280. For the ILM we used the same model studied 
in |l| including the standard density of N/V = lfrn . 
Additional details can be found in [llj . 

We first check that the result © , derived for the non- 
chiral RBM, agrees with the chiral result in the bulk. 
The analytical formula is actually divergent for a = 3/4 
and needs a cutoff in the spectral determinant D(s,g). 
If we calculate the number variance which is given by 
S 2 (L) = L+2 (L — s)R2(s)ds we find that it is strongly 
cutoff dependent. Instead we look at the spectral rigidity 
A 3 (L) = ^r/ L (i 3 - 2l?x + x 3 )Z 2 (x)dx which is not 
sensitive to the cutoff. 

In Figure ^ we show the spectral rigidity in the bulk 
of the chiral RBM obtained from numerical simulation 
at b = 1 along with the analytical formula of the non- 
chiral model for different values of g. The values of g 
were chosen by eye to match the numerical results for 
the corresponding N. We find good agreement with the 
theoretical expectation g 1.17a//V within about 2%. 

In Figure [3 we compare the spectral rigidity in the 
bulk of the ILM with analytical predictions of the RBM. 
Again the values of g for each volume were chosen by eye 
to closely match the ILM. The analytical results using 
these values of g agree well with the ILM up to a scale of 
about 30% of the number of positive eigenvalues (N/2). 
The values of g scale closely to g sw 1.23>/N except for 
the smallest volume (N = 320) where a 5% deviation 
is observed. This corresponds to b w 1.1. The observed 
scaling of the dimensionless Thouless energy in the chiral 



4 



• ILM N= 1280 

RBM N = 1280 b = 0.55 
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FIG. 3: Number variance Y?(L) close to the origin. Points 
correspond to numerical simulations of the ILM. Lines cor- 
respond to numerical simulations of the chiral RBM for the 
given values of b. 

RBM is in agreement with theoretical predictions 0, EH ] 
and lattice results for QCD. 

We do not yet have analytical results for the spectral 
correlations of the RBM at the origin and therefore rely 
solely on numerical simulations. In Figure [3] we show the 
number variance near the origin for the ILM and chiral 
RBM. The corresponding values of b for the RBM were 
also chosen by eye to provide a good fit. The agree- 
ment between both models is very good up to about 50 
eigenvalues but we could not fit all volumes with a sin- 
gle b. The reason is that the chiral RBM has the same 
g ~ >/~N scaling as in the bulk while the instanton liquid 
shows a weaker volume dependence. We find the scaling 
b w 1.4 7V" - 13 which gives g « 1.4 7V - 43 for the ILM. It 
would be interesting to compare these results with lat- 
tice simulations where a g ~ y/~N scaling has also been 
reported close to the origin. 

Finally we mention how the chiral RBM should be 
modified to describe QCD at finite temperature. As usual 
in field theory, temperature is introduced by compacti- 
fying one of the spatial dimensions. Thus the effect of 
temperature in Euclidean QCD is to reduce the effective 
dimensionality of the system to three. Now since the 
effective dimension of the space matches the power-law 
decay of the QCD Dirac operator (~ 1/-R 3 ) one expects, 
according to |J @ , multifractal wavefunctions typical of 
a metal-insulator transition. The same chiral RBM pro- 
posed in this paper may be used in this situation but 
with a = 1 (see [i^ for a model with similar spectral 
properties). The above arguments suggest that if the 
restoration of chiral symmetry at finite temperature is 
dominated by instantons, the physical mechanism lead- 
ing to the quark-gluon plasma state of matter would be 
similar to a metal-insulator transition. Clearly further 
work is needed to explore this exciting relation. 

To conclude, we have proposed a chiral random banded 
model with power-law decay in order to describe the spec- 



tral correlations of the QCD Dirac operator beyond the 
Thouless energy. We have thus combined the asymptotic 
power-law tail observed in instanton liquid models with 
the random matrix approach valid for small spacings. We 
have provided numerical evidence that the resulting chi- 
ral RBM does (at least for the two-point function) de- 
scribe the spectral correlations of the QCD Dirac oper- 
ator well beyond the Thouless energy. Finally we have 
mentioned that at the finite temperature chiral restora- 
tion transition the appropriately modified chiral RBM 
predicts a metal-insulator behavior including multifractal 
wavefunctions and the physics of the Anderson transition 
in the QCD vacuum. 
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sions. A.M.G. was supported by the EU network "Math- 
ematical aspects of quantum chaos". J. CO. was sup- 
ported in part by NSF PHY 01-39929. 
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